Autocorrelation functions in 3D Fully Frustrated Systems 
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D ■ Abstract 

We present a numerical study of autocorrelation functions of a 3D Fully 
Frustrated Ising model (FFIM) simulated by spin-flip Monte Carlo dynamics 

a 

finding simple exponential decay for all the temperature above the critical 
temperature T c for the autocorrelation of squared magnetization and stretched 
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exponential decay for the energy autocorrelation below a temperature T* 
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with T c < T* < T p where T p is the Kasteleyn-Fortuin and Coniglio-Klein 
percolation temperature. The results are compared to those on 2D FFIM to 



> 

On 
O 

in 

[~^. . make light on the relevant mechanism in the onset of stretched exponential 
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relaxation functions. 
Autocorrelation functions of the form 



X3; /(t) = /oexp-(t/r)" (1) 

o 

called stretched exponential or Kohlrausc- Williams- Watts (KWW) function are observed in 
glassy systems and spin glass (SG) well above the respective transition temperature and the 
origin of this kind of behaviour is still matter of debate. 

For example in SG Randeria et al. ]I|] suggested that the onset of the non exponential 
relaxation function coincides with the critical temperature of the ferromagnetic Ising model 
and base their conjecture on the presence of non frustrated ferromagnetic- type clusters of 
interactions. 



Furthermore Huse and Fisher || on the basis of droplet model and more recently Stauffer 
|3|] on the basis of percolation arguments have predicted a KWW autocorrelation function 
also for ferromagnetic Ising model (FIM) in 2D for T <T C where T c is the critical tempera- 
ture. These arguments show that no such behaviour is expected in more than two dimensions 
and Stauffer has corroborated his prediction with simulations in 2 and 3D. 

It is useful to note at this point that in FIM the critical temperature T c coincides with the 
percolation temperature T p of the Kasteleyn-Fortuin and Coniglio-Klein (KF-CK) clusters. 

i 

These clusters can be obtained introducing bond with probability p — 1 — e ~ 2J /( k BT) ^ e _ 
tween nearest-neighbor pair of spins satisfying the interaction (J is the interaction strength, 
ks the Boltzmann constant and T the temperature) and represent correlated spins regions. 

It would be interesting to investigate whether such prediction also holds for frustrated 
systems where the percolation temperature and the critical temperature do not coincide. 

In particular in this work we considered the fully frustrated Ising model (FFIM). Here 
we resume the numerical results on 2D || and present those in 3D. [|5| 

The FFIM is a model with no disorder (unlike the SG and like the FIM) and with the 
interactions assigned in a regular way to give rise to frustration, i.e. a competition between 
them (like the SG and unlike the FIM). [g] 

The FFIM Hamiltonian is given by 

H = -JJ2(^S 1 S 1 -1), (2) 

where e^ = ±1 are quenched regular assigned variables that give rice to frustration and 
Si = ±1 are the Ising spins. || In this model the KF-CK clusters are such that they do 
not include a frustrated loop, that is a closed path of bonds which contains an odd number 
of interactions with e^ = — 1. These clusters overlap regions of not frustrated interactions 
and their percolation temperature is greater then the critical temperature. In Ref. || it is 
studied a general Hamiltonian, with Ising spins and s-state Potts spins, that has the FFIM 
as a particular case for s = 1. It is possible to see that for all s =£ 1 T p of KF-CK clusters 



corresponds to a thermodynamic transition of a s-state Potts model and that this thermal 
transition disappears for s — 1, i.e. for FFMI. This suggests that the KF-CK clusters may 
have an important role also for FFIM. 

We simulate this model using the standard spin flip [|9] dynamics. 

In the 2D FFIM where T p > T c = recent results || shows that stretched exponential au- 
tocorrelation functions are observed for T < T p . Therefore the same qualitatively behaviour 
of 2D FIM is observed but in the FFIM is clear that the onset of stretched exponentials 
coincides, within numerical precision, with T p ^ T c , while in FIM is T p = T c . 

To put more light on the problem we have done simulations on 3D FFIM. 

To this goal we have calculated the autocorrelation function of squared magnetization 
M 2 and of energy E, where for a generic observable A the relaxation function is defined as 

(SA(t)5A{0)) 
fA(t) = {{6A)2) , (3) 

where 5A(t) = A(t) — (A), t is the time and (■) is the thermal average. 
We fit our data with a general function 

m = ft !OL±£Z (4) 

where /o, r, (3 and x are parameters. 

In Fig.|l| we summarize the results for systems with linear size L = 20 and for tempera- 



tures below the estimated |1(J T p (L) of KF-CK clusters in cubic FFIM. 

For f M 2 the estimated parameters are /3 ~ 1 and x ~ for all simulated temperatures 
below T p (L). For f E there is a temperature T*, between T p (L) and T c ~ 1.355 J/ k B (the 
estimated critical temperature [pT|), below which is (3 < 1 and x weakly changes. Simulations 
on L = 30 -T- 50 confirm these data. 

Therefore the 3D FFIM results show a simple exponential decay in M 2 correlation 
functions but a stretched exponential decay in E correlation functions for y < y* with 
T c < T* < T p . In conclusion the results on f E in 3D FFIM, even if less clear than those 
on 2D FFIM, confirm the existence of an onset T* of stretched exponential for the relax- 
ation function well above the critical temperature T c and almost equal (or slightly below) to 



the KF-CK percolation temperature T p that, in a generalized model with also s-state Potts 
variables with s^ 1, corresponds to a thermodynamic transition of the Potts variables. M 
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FIG. 1. Correlation functions for L = 20: (a) j M i and (b) Je at (from left to right) 
ksT I J = 3.5, 2.5, 1.9, 1.8 from respectively 2,7,8,7 independent runs with statistics of 5 • 10 6 
L 3 equilibrium update after 5 • 10 3 L 3 update of thermalization; insets: the corresponding esti- 
mated (3 and x in eq.(4) vs. T/T P (L) where k B T p (L = 20)/ J ~ 3.95. [10] 



